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Results on A, A, H~, and H + production in central Pb+Pb reactions at 2CL4, 30^4, 40A, 80A, 
and 158A GeV are presented. The energy dependence of transverse mass spectra, rapidity spectra, 
and multiplicities is discussed. Comparisons to string hadronic models (UrQMD and HSD) and 
statistical hadron gas models are shown. While the latter provide a reasonable description of all 
particle yields, the first class of models fails to match the H~ and H + multiplicities. 

I. INTRODUCTION be formed. In this state of matter the normal confine- 

ment of quarks and gluons in hadrons is removed and the 
partons can exist as quasi-free particles in an extended 
Heavy ion reactions at ultra-relativistic energies allow region of space-time. One of the first signatures proposed 
the study of strongly interacting matter at extreme tern- for the formation of a QGP state was an enhancement of 
peratures and densities. It is expected that under such strange particle production in A+A with respect to p+p 
conditions eventually a quark gluon plasma (QGP) will collisions [1]. The argumentation relies on the assump- 
tion that gluon fusion processes, which may be dominant 

in a QGP, produce additional ss pairs[2|. This in turn 

will cause the abundance of strange quarks to reach its 
'deceased chemical equilibrium value in much shorter times than 



would be possible in a pure hadronic scenario. In fact, 
the enhanced production of strange particles has been ob- 
served already quite early in high energy nucleus-nucleus 
collisions [3, lj| • It has also been demonstrated that the 
enhancement is most pronounced for the multiply strange 
hyperons 5 and fl 0, @, 0| ■ 

However, systematic studies of hadron production in 
nucleus-nucleus collisions have shown that strangeness 
enhancement is not only seen at high energies, such as 
top SPS and RHIC energies, but it is also observed at 
lower energies (^/s NN < 5 GeV) Q where no QGP forma- 
tion is expected. Actually, here the production of A and 
5 exhibits an even stronger enhancement than present 
at top SPS or RHIC energies [KMO- Generally, it is 
found that the abundances of strange particles in central 
A+A reactions are similar to those expected from sta- 
tistical hadro n g as models assuming a grand canonical 
ensemble [l2|, [l3[ . While the enhancement at lower ener- 
gies can to a certain extent also be explained by trans- 
port models, at higher energies additional mechanisms 
have to be involved in order to reach chemical equilib- 
rium values via a dynamical evolution. It has, e.g., been 
suggested that multi-pion reaction processes can lead to 
an accelerated equilibration of anti-hyperon production 
in nucleus-nucleus collisions [14]]. Especially at larger 
densities, as present close to the QGP phase boundary, 
processes like this might drive the system quite fast to a 
chemical equilibrium state [l||. Still it is an open ques- 
tion whether such dynamical explanations are applicable 
as well at lower energies. On the other hand, it was 
suggested that particle production via strong interaction 
always follows the maximum entropy principle and there- 
fore hadron abundances are naturally close to the out- 
come of statistical processes [la, [ill [01 LLSl • The mea- 
surement of hyperon production in an energy range below 
top SPS energy (^/s NN < 17.3 GeV) provides important 
constraints on both, the statistical and transport model 
approach. Recent results on kaon production in central 
nucleus-nucleus collisions [201 ] indicate a sharp maximum 
of the ratio (K + )/(ir + ) and a sudden change in the energy 
dependence of the (m t ) — rriQ of pions, kaons, and protons 
at a beam energy of 30>1 GeV, where m t = \Jvl + TO o ^ s 
the transverse mass, nig the rest mass and pt the trans- 
verse momentum. These anomalies can be interpreted 
as a signal for the onset of deconfinement [2l|, [22| and 
might also be visible in the energy dependence of hy- 
peron yields. 

The data discussed here represent an extension of pre- 
viously published results [10, [ill, (23| in order to provide a 
complete study of the energy dependence of hyperon pro- 
duction at the CERN-SPS. Some of the data discussed 
here have been presented as preliminary before [7L [24 [2a] . 
However, the data shown in this publication are the re- 
sult of a completely new and independent analysis, which 
treats all datasets in a consistent manner. In particular, 
the results for A and A include a correction for the feed- 
down from weak decays, which was not applied in the 
previous publication [101 ] . 



TABLE I: Summary of the analyzed datasets. The central- 
ity fraction corresponds to the most central part of the total 
inelastic cross section. The Glauber model was used to de- 
termine the averaged number of wounded nucleons per event 

(N w ). 
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II. DATA ANALYSIS 
A. Experimental setup and data sets 

The data were taken with the NA49 large acceptance 
hadron spectrometer at the CERN SPS. A detailed de- 
scription of the apparatus can be found in [2g]. With 
this detector, tracking is performed by four large-volume 
Time Projection Chambers (TPCs) in a wide range of 
phase space. Two of these are positioned inside two su- 
perconducting dipole magnets. In order to assure a simi- 
lar detector acceptance for all datasets, the magnetic field 
was chosen proportional to the beam energy. A measure- 
ment of the specific energy loss dE/dx in the TPC gas 
with a typical resolution of 4% provides particle iden- 
tification at forward rapidities. Time-of-flight detectors 
improve the particle identification at mid-rapidity. Cen- 
tral Pb+Pb reactions were selected by imposing an upper 
threshold on the energy measured in the projectile frag- 
mentation region. For this measurement the Zero Degree 
Calorimeter (ZDC) was used which is positioned down- 
stream of the TPCs. A collimator in front of the ZDC 
assures that the acceptance of the calorimeter matches 
the phase space of the projectile fragments and specta- 
tor nucleons. 

We present in this paper an analysis of central Pb+Pb 
events taken at beam energies of 2QA, 30A, 40A, 80A, 
and 158A GeV in the years between 1999 - 2002. The 
properties of the different datasets are summarized in 
Tabled] The 158^4 GeV dataset has an online centrality 
trigger on the 23.5% most central events, of which the 
10% most central were selected offline. 



B. A (A) and H (H + ) reconstruction 

A and A hyperons were found by reconstructing their 
charged deca ys A — > ir~ + p and A — > tt + + p (branching 
ratio 63.9 % [23]). In a first step pairs were formed of all 
positively charged particles with all negatively charged 
ones. Their tracks were reconstructed by a global track- 
ing algorithm that connects the track parts registered 



xitr 



x10 d 



x10 J 




1.1 1.15 1.2 

(PX) (GeV/c 2 ) 
x10 3 



20A GeV 

(e) 
N 



1.1 1.15 1.2 

(p.rc) (GeV/c 2 ) 
x10 6 



158AGeV 
(0 



0.01 




1.1 1.15 1.2 

(p,it + ) (GeV/c 2 ) 

x10 3 



£0.4 



0.3 



0.2 



0.1 



u 1.3 1.4 
(A,ji") (GeV/c 2 ) 



1.3 1.4 

(A,ji") (GeV/c 2 ) 



30A 

'(g) 


GeV 

__+ 


aJ 


V 


* 


H 



1.1 1.15 1.2 

(p,jt + ) (GeV/c 2 ) 
x10 3 



1 58A GeV 

4(h) 



1.3 1.4 

(A\ji + ) (GeV/c 2 ; 




1.3 1.4 

(A",ji + ) (GeV/c 2 ) 



FIG. 1: The invariant mass distributions of all A, A (upper 
row), E~, and H + (lower row) candidates in central Pb+Pb 
collisions at the lowest and highest analyzed beam energies. 
The full curves represent a fit to signal and background as 
described in the text, while the dashed curves show the back- 
ground contribution only. The gray vertical lines denote the 
PDG masses [H. 
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FIG. 2: The total reconstruction efficiency of A (A) (upper 
row) and H~ (H + ) (lower row) as a function of pt (left pan- 
els), and as a function of rapidity (right panels) for central 
Pb+Pb at 158A GeV. The symbols denote the efficiency for 
the standard analysis procedure (STD). In addition, the A 
(H~) efficiencies for two other selection criteria are shown 
(dashed: CUT-A, dotted: CUT-B, see text). 



in the different TPCs. Only tracks with more than 10 
reconstructed points were accepted. By requiring a dis- 
tance of closest approach (DCA) between their trajecto- 
ries of less than 0.5 cm anywhere between the position 
of the first measured points on the tracks and the tar- 
get plane, V° candidates were identified. A set of ad- 
ditional cuts was imposed in order to reduce the com- 
binatorial background due to uncorrelated pairs. Iden- 
tification of (anti-)protons via their specific energy loss 
(dE/dx) in the TPCs reduces the contribution of pairs 
with a wrong mass assignment. The measured dE/dx 
was required to be within 3.5 standard deviations from 
the predicted Bethe-Bloch value. A A (A) candidate was 
accepted if the reconstructed position of its decay ver- 
tex is at least 30 cm downstream from the target and 
outside the sensitive volume of the TPCs, to avoid inef- 
ficiencies resulting from an insufficient separation of the 
clusters of the two tracks. The trajectories of the A (A) 
candidates were extrapolated back to the target plane 
to determine their impact parameters 6 X (in the mag- 
netic bending plane) and b y relative to the interaction 
point. Non-vertex candidates were rejected by requir- 
ing |6 X | < 0.5 cm and |6 y | < 0.25 cm. A further re- 
duction of the combinatorial background was achieved 
by placing a requirement on the angle 9* between the 
flight direction of the mother particle and of its positive 
daughter, determined in the center-of-mass system of the 
A (A) candidate: -0.95 < cos6>* < 0.75. Figure up- 
per row, shows as examples the resulting invariant mass 



spectra at 20A and 158,4 GeV for A and A. The posi- 
tion of the peaks in the m; n v distribution agrees with the 
nominal A mass determined by the particle data group 
|27| . From a fit with a Gaussian typical mass resolutions 
of er m « 2 MeV/c 2 are determined, which depend only 
slightly on phase space and beam energy. Generally, the 
signal to background ratio (S/B) is worse for A than for 
A, due to the lower yield of real A. While S/B decreases 
with energy for A, it is constant for A. However, the 
shape of the combinatorial background depends on beam 
energy in both cases. 

S~ (S + ) candidates were identified via the decay chan- 
nel S~ — > A + 7r~ (S + — > A + 7r + ) which has a branching 
ratio of 99.9 % [27J]. To reconstruct the S~ (E+), A 
(A) candidates were selected in an invariant mass win- 
dow of 1.101 - 1.131 GeV/c 2 and combined with all 
measured negatively (positively) charged particles in the 
event. The A (A) candidates were required to pass the 
same cuts as described above, with the exception of the 
cuts on b x / y and cos 6* , which were not applied here. The 
reconstructed S~ (S + ) candidates should point back to 
the interaction vertex, while the pions from A (A) and 
the S~ (2 + ) decay will on average have a larger im- 
pact parameter. To reject non-vertex candidates, upper 
limits of |6 X | < 0.5 cm and |6 y | < 0.25 cm were there- 
fore imposed on the S~ (S + ) candidates. Pions com- 
ing from the primary interaction point were removed by 
a cut of | fry | > 1.0 cm for the negatively (positively) 
charged tracks associated to the S~ (S + ) decay ver- 



tex and |6 y | > 0.5 cm for the negatively (positively) 
charged daughter tracks of the A (A) candidates. An 
additional improvement of the signal to background ra- 
tio was achieved by requiring that the A decay vertex 
and the pion track were measured on the same side of 
the TPCs relative to the beam pipe. The lowest beam 
energy where a significant S + signal could be extracted 
is 30A GeV, while S~ could be analyzed at all available 
energies. Figure [TJ lower row, shows the invariant mass 
distributions for S~ and S + candidates after all selec- 
tion criteria at the lowest and highest available energies, 
respectively. Similarly as in the case of the A (A) an 
excellen t ag reement of the peak positions with the PDG 
masses [27j is observed. The typical mass resolution, as 
obtained from a fit with a Gaussian is a m « 4 MeV/c 2 . 
The dependence of the shape of the combinatorial back- 
ground on the beam energy is less pronounced than in 
the case of A (A) . 

The invariant mass spectra were fitted to the sum 
of a polynomial and a signal distribution, determined 
from the simulation procedure described below. The raw 
yields of A, A, H~, and S + are obtained by subtract- 
ing the fitted background and integrating the remaining 
signal distributions in a mass window of ±11 MeV/c 2 
(±10 MeV/c 2 ) around the nominal A (S) mass. 



C. Correction for acceptance and reconstruction 
inefficiency 

Detailed simulations were made to correct the yields 
for geometrical acceptance and losses in the reconstruc- 
tion. For this purpose, samples of A and S were generated 
in the full phase space accessible to the experiment with 
rat spectra according to: 
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and Gaussian shaped distributions in rapidity y. The 
Geant 3.21 package [28| was used to track the generated 
particles and their decay products through the NA49 de- 
tector. Dedicated NA49 software was used to simulate 
the TPC response taking into account all known detec- 
tor effects. The simulated signals were added to those 
of real events on the raw data level and subjected to the 
same reconstruction procedure as the experimental data. 
The acceptances and efficiencies were calculated in bins 
of pt (mt — rao) and y as the fraction of the generated A 
(S) which traverse the detector, survive the reconstruc- 
tion and pass the analysis cuts. Of all produced hyperons 
~ 50% (40%) of the A (5) appear in the acceptance of the 
detector, i.e. all decay particles are seen in the sensitive 
detector volume. The reconstruction algorithm together 
with the cuts to suppress the combinatorial background 
reduce this fraction further to ~ 6% (4%) at 158A GeV. 
In addition, inefficiencies due to the high track multiplic- 
ity cause a further reduction. At 158A GeV this effect 
is most pronounced and reduces the integrated efficiency 
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FIG. 3: The feed-down contribution to A (A) around mid- 
rapidity as a function of pt (left panel), and as a function 
of rapidity (right panel) for central Pb+Pb at 158A GeV. 
The symbols denote the feed-down for the standard analysis 
procedure (STD). In addition, the feed-down to A for two 
other selection criteria are shown (dashed: CUT-A, dotted: 
CUT-B, see text). 



to ss 2% (1%) for A (S). At lower energies the influence 
of the occupancy is weaker and thus the total efficiency 
increases to sw 4% (2%) at 20,4 GeV. The upper row of 
Fig. [2] shows the total reconstruction efficiency which in- 
cludes acceptance and all reconstruction inefficiencies for 
A and A at the highest beam energy, where the effects of 
the high track density are largest. Also included in Fig. [2] 
are efficiencies that have been calculated for two analy- 
sis strategies different from the default version described 
above. The first one (CUT-A, shown as dashed lines) 
employs a set of selection criteria that depend on the 
sub-detector in which a V was measured and which were 
optimized for a large signal-to-background ratio [23| ■ The 
second strategy (CUT-B, shown as dotted lines) uses the 
same cuts as described above, but in addition only ac- 
cepts tracks which lie outside the high track density re- 
gion (i.e. 4 cm above or below the middle plane of the 
TPCs). This criterion allows to minimize the losses due 
to the high occupancy at the expense of a much reduced 
acceptance in particular at low pt- It was used in a pre- 
vious analysis of the A (A) at 158 A GeV published in 
[lCj . Both approaches drastically reduce the number of 
reconstructed A (S). Therefore, they were not used as 
the standard procedure in this analysis, but can serve 
as a cross check that helps to estimate systematic errors 
(see section Hi Ej) . 



D. 



Correction of feed-down to A (A) 



The measured yield of A and A contains, in addition 
to the directly produced particles, contributions from the 
decay of heavier hyperons. The A (A) resulting from 
electromagnetic decays of E° (E°) cannot be separated 
from the directly produced ones. Thus the here presented 
yields always represent the sum A+E° (A+E°). The con- 
tribution to A (A) from weak decays, however, depends 
on the chosen analysis cuts, since these decay products 
originate from decay vertices with a sizable distance from 



TABLE II: Summary of the systematic errors on the 
dN/dy values. 

Background Efficiency pt Extra- Feed-down Quadratic 
subtraction correction polation correction sum 
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the main interaction point. Since the NA49 acceptance 
for A (A) favours those that decay at large distances, the 
contribution of feed down A (A) can be quite substantial. 
Therefore, we have calculated a correction for the feed- 
down from S~ + S° (S+ + S°) decays to the measured 
A (A) sample with the same simulation procedure as de- 
scribed above for the efficiency correction. In this case a 
sample of S~ and S° (S + and S°) was generated as input 
to the NA49 simulation chain. The feed-down correction 
was then calculated in bins of pt (wit — tio) and y as 
the fraction of reconstructed A (A) which originate from 
E~ + El° (S + + S°) decays and pass the same analysis 
cuts. The input distributions and yields of the S~ (S + ) 
are the ones measured by NA49 and presented in this 
publication. For the S° (S°), which are not measured, 
the same phase space distributions were assumed. Their 
yields are calculated from the ones of S~ (S + ) which are 
scaled by the S°/S~ (S°/S + ) ratios taken from a statis- 
tical model fit [29|] . Figure [3] shows as an example the 
calculated feed-down contribution to A (A) as a function 
of pt and rapidity. The feed-down is largest at low p t and 
mid-rapidity and larger for A (20 - 30 % at 158A GeV) 
than for A (5 - 15 % at 158A GeV). While for A no 
significant dependence of the feed-down on the beam en- 
ergy is observed, the feed-down to A reduces to 3 - 8 % at 
20 A GeV. Also included in Fig. [3] are the feed-down con- 
tributions to A for the two alternative analysis strategies 
described in the previous section (dashed line: CUT-A, 
dotted line: CUT-B). Since the fraction of S seen in the 
reconstructed A sample depends on the selected analysis 
cuts, the feed-down contribution has to be evaluated for 
each approach separately. 



E. Systematic errors 

There are several contributions to the systematic error 
which are summarized in Tabic [Til One of them results 
from uncertainties in the determination of the combina- 
torial background. This uncertainty can be estimated by 
varying the degree of the polynomial used to fit the back- 
ground and the invariant mass range in which the fit is 
performed. It is found that this systematic error is 3% 
for A and 5. 

Another contribution arises from imperfections in the 
description of the detector response by the simulation 
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FIG. 4: The differences between the fully corrected results 
of the standard procedure and of the two alternative anal- 
ysis strategies (see section III Ef) for A (upper row) and H~ 
(lower row) in central Pb+Pb at 158A GeV. Shown are the 
pt dependence at mid-rapidity (left panels) and the rapid- 
ity dependence (right panels). The gray boxes illustrate the 
systematic error estimate. 



procedure which result in systematic uncertainties in the 
efficiency calculation. It was verified that all distribu- 
tions of geometrical and kinematical parameters that are 
relevant in the reconstruction procedure (see section. Hl Bp 
are in a gre ement between simulated and measured data 
|30t l3ll . l32j . Still there can be remaining discrepancies 
which constitute a source of systematic error. Its mag- 
nitude can be estimated by varying the selection criteria 
in the analysis procedure and checking the consistency of 
the final result. This was done, e.g., by comparing data 
points obtained with the alternative analysis strategies 
described in section [ETC] (CUT-A and CUT-B) to the re- 
sults for the standard analysis (see Fig.@}. Shown are the 
differences AJV = V(STD)- iV(CUT-A(B)) as a function 
of pt and rapidity, both for A and S - . Even though the 
efficiencies are lower by almost a factor 2 in some regions 
of phase space (see Fig. [5]) and are subject to different 
systematic effects (e.g. influence of high track density or 
background) the results are in agreement. The deviations 
are consistent with a systematic error of « 10% for A and 
E~~ at all beam energies (see Fig. 2]). Additionally to the 
studies presented in Fig. [4j a further investigation was 
performed in order to test whether the A reconstruction 
is sensitive to the cut applied to the DCA. For this pur- 
pose the DCA-cut was relaxed to 1.5 cm (default: 0.5 cm) 
and the result of this analysis was compared to the stan- 
dard procedure. It was found that the deviations between 
the two approaches also agree with the systematic error 
estimate given in Table [Til 



TABLE III: The rapidity densities dN/dy at mid-rapidity (A/A: \y\ < 0.4, H"/H + : \y\ < 0.5), the total multiplicities (TV), 
the RMS-widths of the rapidity distributions RMSy calculated from the fits shown in Fig. [9] the averaged transverse masses 
(rat) — mo, and the inverse slope parameters T at the different beam energies E^ cam . The first error is statistical, the second 
systematic. 





-^bcam 


Cent. 


dN/dy 


(N) 


RMSy 


(m t ) - m 


T 




(A GeV) 


(%) 








(MeV/c 2 ) 


(MeV) 


A 


20 


7 


13.4±0.1±1.1 


27.1±0.2±2.2 


0.70±0.01±0.06 


297± 4±24 


244± 3±12 




30 


7 


14.7±0.2±1.2 


36.9±0.3±3.3 


0.89±0.02±0.08 


310± 5±25 


249± 2±13 




40 


7 


14.6±0.2±1.2 


43.1±0.4±4.3 


1.11±0.08±0.11 


327± 5±27 


258± 3±13 




80 


7 


12.9±0.2±1.0 


50.1±0.6±5.5 


1.28±0.02±0.14 


338± 7±27 


265± 4±13 




158 


10 


9.5±0.1±1.0 


44.9±0.6±8.0 


— 


368± 7±28 


301± 4±15 


A 


20 


7 


0.10±0.02±0.01 


0.16±0.02±0.03 


0.62±0.14±0.14 


407±72±47 


339±56±31 




30 


7 


0.21±0.02±0.02 


0.39±0.02±0.04 


0.69±0.05±0.08 


357±32±30 


284±13±26 




40 


7 


0.33±0.02±0.03 


0.68±0.03±0.07 


0.77±0.05±0.08 


371±22±31 


301±10±27 




80 


7 


0.82±0.03±0.08 


1.82±0.06±0.19 


0.83±0.05±0.09 


363±19±30 


292±10±27 




158 


10 


1.24±0.03±0.13 


3.07±0.06±0.31 


1.00±0.03±0.09 


388±13±31 


303± 6±27 


3 


20 


7 


0.93±0.13±0.10 


1.50±0.13±0.17 


0.64±0.08±0.07 


289±27±29 


221±14±13 




30 


7 


1.17±0.13±0.13 


2.42±0.19±0.29 


0.73±0.14±0.09 


278±19±28 


233±11±14 




40 


7 


1.15±0.11±0.13 


2.96±0.20±0.36 


0.94±0.13±0.11 


285±17±29 


222± 9±13 




80 


7 


1.22±0.14±0.13 


3.80±0.26±0.61 


0.98±0.25±0.16 


317±22±32 


227±14±14 




158 


10 


1.44±0.10±0.15 


4.04±0.16±0.57 


1.18±0.18±0.17 


327±13±33 


277± 9±17 


H + 


20 
30 


7 

7 


0.05±0.01±0.01 


0.12±0.02±0.03 


0.76±0.35±0.17 


326±60±33 


311±75±31 




40 


7 


0.07±0.01±0.01 


0.13±0.01±0.02 


0.65±0.13±0.09 


337±36±34 


277±32±28 




80 


7 


0.21±0.03±0.02 


0.58±0.06±0.13 


0.87±0.29±0.20 


298±38±30 


255±23±26 




158 


10 


0.31±0.03±0.03 


0.66±0.04±0.08 


0.73±0.08±0.09 


384±26±38 


321±15±32 



In case of A and A also the uncertainties in the feed- 
down contribution have to be taken into account. Here, 
the errors of the measurements of spectra and yields of 
S~ and S + translate into a systematic error caused by 
the feed-down correction. For A this error is small (1%), 
since the correction itself is not too substantial and the 
S~ measurement is relatively accurate. In case of A, 
however, the larger feed-down contribution and the larger 
errors of the S + data also result in a larger systematic 
error of 7%. 

In the range of the errors the data presented here agree 
well with the previously published results where available 
[lOL [23j | . The differences compared to the A and A yields 
given in [101 ] are due to the feed-down contribution which 
has not been subtracted from the old results. 

Since for the A (A) the range down to p t = GeV/c is 
measured in most of the rapidity bins, the systematic ef- 
fects due to extrapolations into unmeasured p t regions is 
negligible. Only in those y bins where extrapolations are 
necessary an additional systematic error of 4% is added 
in quadrature. However, for the S~ (S + ) analysis, this 
is introducing an additional systematic error in the full 
range of the rapidity distributions. It was estimated by 
using different assumptions for the spectral shape. The 
standard approach (fit to an exponential, see section HV)) 
was compared to a fit with a hydrodynamically inspired 



blast wave model [32 
found to be 3%. 



The difference on the dN/dy was 
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FIG. 5: The transverse mass spectra of A (left panel) and A (right panel) at mid-rapidity (\y\ < 0.4) for 5 different beam 
energies. The data points are scaled for clarity. Only statistical errors are shown. The solid/dashed lines represent a fit with 
an exponential, where the solid part denotes the mt range in which the fit was performed. The dotted lines are the results of 
a fit with a blast wave model [33J ] (see text for details). 
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FIG. 6: The transverse mass spectra of H~ (left panel) and H + (right panel) at mid-rapidity (\y\ < 0.5) for 5 (4) different beam 
energies. The data points are scaled for clarity. Only statistical errors are shown. The solid/dashed lines represent a fit with 
an exponential, where the solid part denotes the mt range in which the fit was performed. The dotted lines are the results of 
a fit with a blast wave model [33j ] (see text for details). 
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The extraction of the total multiplicities requires in 
addition an extrapolation into the unmeasured rapidity 
regions. The systematic error that is introduced by this 
extrapolation depends on the beam energy, since the frac- 
tions of the longitudinal phase space covered by the mea- 
surements also change with energy. Also, the shape of the 
y spectra is not always very well determined, especially 
for A at 80A and 158A GeV. By using different assump- 
tions for the spectral shape in the unmeasured region, 
as defined in Eq. fl} and Eq. @, the additional sys- 
tematic error on the total multiplicities was estimated. 
For A a variation between 1% at 20A GeV and 14% at 
158 A GeV was obtained and for A this systematic error 
is largest at the lowest energy (20%) and decreases to 2% 
at 158A GeV. In case of the S~ this contribution ranges 
between 2% at 20A GeV and 12% at 80A GeV, while for 
the 5+ it is between 5% (158 A GeV) and 20% (30A and 
80 A GeV). 



III. TRANSVERSE MASS SPECTRA 
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FIG. 7: The (m t )-mo values for central Pb+Pb and Au+Au 
reactions as a function of ,/s NN . The systematic errors are 
represented by the gray boxes. Filled symbols correspond to A 
and H~ , while open symbols denote A and H + . Also shown are 
data from the NA57 collaboration [H3, from AGS [HHI 
and RHIC experiments [371 . l3a . |391 |. as well as p+p data on 
A [4011. The lines are calculations with the UrQMDI. 3 model 



The transverse mass spectra of A and A measured 



around mid-rapidity (\y\ < 0.4) are shown in Fig. [5] and 
the ones of S~ and S + (\y\ < 0.5) in Fig.[6j The mt spec- 
tra were fitted by an exponential in the transverse mass 
range m t — m > 0.2 GeV/c 2 as defined in Eq. ([T]). The 
resulting inverse slope parameters T are summarized in 
Table UTTl Due to the significant curvature of the mt spec- 
tra equation Eq. ([T|) does not provide a satisfactory de- 
scription of the data over the whole m t range. Therefore, 
the rrit spectra were additionally fitted by a blast wave 
model which assumes a transversely expanding emission 
source [33j . The parameters of this model are the freeze- 
out temperature Tf and the transverse flow velocity f3 a 
at the surface. Assuming a linear radial velocity profile 
f3t(r) = p s r/R, which is motivated by hydrodynamical 
calculations, the m t spectrum can be computed from 

1 dJV f R (p t smhp\ fm t cosh p 
— a r K rdrm t l [ Kd — 

' r (2) 
where R is the radius of the source and p = tanh pt 

is the boost angle. Since the measurements for the dif- 
ferent particle species do not provide an equally good 
constraint on the fit procedure if both parameters are 
allowed to vary freely, the transverse flow velocity was 
fixed to (/3 t ) = 2/3 f3 s = 0.4. The results of the fits are 
shown as dotted lines in Figs. [5] and [6] and the obtained 
fit parameters Tf are listed in Table ITVl They turn out to 
be significantly lower for S~ than for A at all beam ener- 
gies. This difference is also visible for the anti-particles, 
although less pronounced. Even though this observation 
is based on a relatively simple model, it might indicate 
that the transverse mass spectra of A and S are not de- 
termined by the same kinetic freeze-out condition. 



TABLE IV: The parameter Tf resulting from the fit with the 
blast wave model. (/3t) was fixed in all cases to 0.4. Tf is 
given in MeV. Errors are statistical only. 



Beam energy 


T f ( A) 


Tf(A) 


Tf(H-) 


T f (H+) 


20A GeV 


100T2 


166T38 


82±7 


— 


30A GeV 


107T1 


134T9 


83±5 


122T30 


40A GeV 


115T2 


143T7 


82±4 


127T17 


80A GeV 


121T2 


136T6 


95±8 


108T12 


158A GeV 


140T2 


146T3 


109T5 


156±9 



To allow for a model independent study of the energy 
dependence of mt spectra, the averaged transverse mass 
(mt) — mo was calculated. Since for A, A, and S - es- 
sentially the whole range down to mt — mo = is cov- 
ered, (mt) — mo can be extracted from the data alone. 
However, in order to extrapolate up to a common upper 
limit in mt — mo, fit functions were used as well. For 
this purpose two different fits were used: The blast wave 
model, as shown in Fig. [5l and a fit with a double expo- 
nential (not shown) that also provides a good description 
of the data. The different approaches allow to estimate 
the systematic error. For S + also an extrapolation to 



fn% — mo = is needed. 
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FIG. 8: The transverse mass spectra of A, E _ , A, and H + for central Pb+Pb collisions in different rapidity bins at 20A, 30A, 
40A, 80A, and 158^4 GeV. Every rapidity bin has a width of 0.4. Bin yO starts at -2.0. The data points are scaled for clarity. 
Only statistical errors are shown. 
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FIG. 9: The rapidity spectra of A, E~, A, and H + for 5 different beam energies. The open symbols show data points reflected 
around mid-rapidity. The systematic errors are represented by the gray boxes. Solid lines are fits to the data points, used to 
extract the total yields, while dashed lines are the fits to 2~ and H + from [23(]. The gray area in the A spectrum at 158A GeV 
depicts the uncertainty due to the different extrapolations (see text). 
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The resulting values for (m t ) —Too, corresponding to an 
interval < TOt — mo _ 2 GeV/c 2 , are listed in Table HTT1 

Figure [7] shows the dependence of (m t ) — too on ^/s NN 
for the data presented here in comparison to measure- 
ments done by NA57 at the SPS [& H, b Y E896 and 
E917 at the AGS (Hill and b Y STAR and PHENIX at 
RHIC (13, H|3|. The {m t )-m values derived from the 
NA57 spectra agree with the NA49 results. In the SPS 
energy range only very little variation of (mt) — Too with 
,/s NN is observed, followed by a slight increase towards 
RHIC energies. The (mt) — too of A is generally higher 
by rs 200 MeV/c 2 than the one observed in p+p reac- 
tions [4(| at all center-of-mass energies. For pions, kaons 
and protons a sudden change in the energy dependence 
of (mt) — m around y/s NN =7-8 GeV was observed 
[2Cj. Since currently no data at lower energies are avail- 
able, it cannot be established whether a similar feature 
is present in the energy dependence of (mt) — mo for hy- 
perons. However, the remarkably small energy variation 
shown in Fig. [7] would still be in line with the behavior 
observed for the other particle species. 



TABLE V: The parameter a and j/o resulting from the fits 
with the sum of two Gauss functions (see Eq. (|3])) to the 
rapidity distributions of A and S - . 



Beam energy °"(A) 



Vo(A) 



q(5 ) y (5 ) 



20A GeV 0.51±0.01 0.49T0.01 0.45T0.08 0.45±0.07 

30A GeV 0.66T0.02 0.59T0.01 0.56T0.15 0.47T0.11 

40A GeV 0.91T0.06 0.65T0.04 0.76T0.16 0.54T0.12 

8(L4 GeV 0.87T0.07 0.94T0.06 0.71T0.32 0.68T0.13 

158,4 GeV 1.18±0.18 



The measurements on (m t ) — too are also compared 
to the string hadronic model UrQMD1.3. While this 
model in principle reproduces the observed near indepen- 
dence of (wt) — mo on ^/s NN in the SPS energy region, it 
fails to match its magnitude. The calculation is always 
rs 100 MeV below the data. Additionally, this version 
of UrQMD does not describe the slow increase towards 
RHIC. 



IV. 



RAPIDITY SPECTRA 



Figure [5] summarizes the mt spectra of A, S - , A, and 
S + as measured in different rapidity bins. The data 
points cover a large fraction of the phase space and thus 
allow to extract rapidity distributions by integrating the 
transverse mass spectra. Table IIIII summarizes the re- 
sulting rapidity densities around mid-rapidity and Fig. [§] 
shows the resulting y spectra. For A a clear evolution of 
the spectral shape with beam energy is observed. While 
the rapidity spectrum at 20^4 GeV has an almost Gaus- 
sian shape, a plateau around mid-rapidity is develop- 
ing that widens with increasing energy. At 158A GeV 
the spectrum is finally constant in the measured rapidity 



range. This reflects the continuous change of the rapid- 
ity distribution of the net-baryon number in this energy 
range [43J]. While at lower energies the final state distri- 
bution of the incoming nucleons looks thermal, the ra- 
pidity distribution of the net-baryons develops a distinct 
minimum at mid-rapidity with increasing energy due to 
incomplete stopping. Since A carry a significant frac- 
tion of the net-baryon number they follow this change 
to a large extent. A similar behavior, although less pro- 
nounced, is visible for the S~ as well. A and S + , on the 
other hand, are well described by Gaussians at all beam 
energies. 



00 

rr 




FIG. 10: The RMS widths of the rapidity distributions 
RMSy, normalized by the projectile rapidity j/proj, as a func- 
tion of J/proj- The systematic errors are represented by the 
gray boxes. 

In order to determine total multiplicities, extrapola- 
tions into the unmeasured y regions are needed. There- 
fore, A and S~ were fitted with a sum of two Gauss 
functions placed symmetrically around mid-rapidity 



ON 
dy 



ex exp 



(y - yo) 2 

2a 2 



exp 



(y + yo) 

2a 2 



(3) 



The resulting fit parameters are listed in TablefV] In case 
of the S~ at 158,4 GeV a single Gaussian turned out to 
provide a better fit to the data (solid line in Fig. [9]). For 
the A distribution at 158 A GeV a fit cannot be performed 
since the measurement does not allow to determine the 
end of the dN/dy distribution. Here, the extrapolation 
has to be based on different assumptions on the spec- 
tral shape. An upper limit on the contribution from the 
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unmeasured parts can be derived by using the measured 
net-proton distribution at 158^4 GeV |44j to describe the 
tails. Another approach is to assume the same shape 
for the A rapidity distribution as has been measured for 
central S+S reactions at 200 A GeV [J], which then re- 
sults in a lower total yield. The multiplicity quoted in 
Table [TTT1 is the average between both extrapolations and 
their difference is taken as its systematic error. 
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FIG. 11: The rapidity densities dN/dy at mid-rapidity (A/A: 
\y\ < 0.4, E-/H+: \y\ < 0.5) in central Pb+Pb and Au+Au 
collisions as a function of ^/s NN . The systematic errors are 
represented by gray areas, mostly hidden by the symbols. 
Also shown are data from the NA57 collaboration [g, [3J], 
as well as from AGS [351 13a . l45fl and RHIC experiments 



For A and S+ a single Gauss function was used to de- 
rive the total yields. The resulting fit parameters a are 
identical to the values for RMS y tabulated in Table IIIII 
The rapidity spectra of S~ and S + at 158 A GeV also 
agree well with a fit to the previously published data 
(dashed lines in Fig. [5]). Figure [T01 summarizes the en- 
ergy dependence of the RMS y values. While the widths 
of the A and S + distributions agree with each other and 
exhibit an approximately linear dependence on the pro- 
jectile rapidity j/ pro j (RMS y rs 0.3 j/proj? see dashed line 
in Fig. [TO)), the A and S~ show a different behavior. Here 
RMS y /y pi0i is larger and also clearly energy dependent. 
The effect is more pronounced for the A than for the S - . 



V. PARTICLE YIELDS 

Figure [TTJ shows the rapidity densities around mid- 
rapidity as a function of ,/s NN . The energy dependence 
of dN/dy for A exhibits a complicated structure. It 
rises from AGS to a maximum at a beam energy of 
30A GeV, then drops towards top SPS energy and rises 
again slowly to ^/s^ — 130 GeV. This can be under- 
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FIG. 12: The p/p [43], A/A, H+/3", and «+/0~ llj ratios 
around mid-rapidity (A: \y\ < 0.4, H and Q: \y\ < 0.5) in 
central Pb+Pb and Au+Au collisions as a function of ys NN . 
Also shown as open symbols are data from the SPS experi- 
ments NA44 [13 and NA57 ill, from AGS il [M HE Ell 
and RHIC experiments 0, [M M, [H, [13, H |H2 . The sym- 
bols are slightly displaced for clarity. 



stood by an interplay of the slow rise of the A multiplic- 
ity from .Ebeam = 30A GeV on (see Fig. IT5k ) and the 
pronounced change of shape seen in the rapidity distri- 
bution in the same energy region (see Fig.[H|). Since the 
A yield gets distributed more and more evenly along the 
rapidity axis, the mid-rapidity dN/dy is reduced above 
Ebcum = 30A GeV. At some point the redistribution 
along y is compensated again by the further increase of 
the A multiplicity, so that the rapidity density dN/dy 
is again higher at RHIC. Such a significant structure in 
the energy dependence is not observed for the S - , where 
the mid-rapidity dN/dy increases more smoothly by a 
factor of w 2 from E hcam = 20 A GeV towards RHIC. 
However, also here a small structure in the energy depen- 
dence is visible between 20 A and 80^4 GeV. For A and 
S + , where no change in the shape of the dN/dy spectra 
is seen, the mid-rapidity dN/dy values increase rapidly 
over w 2 orders of magnitude between -Ebcam = 20A GeV 
and y/s^ = 130 GeV. 

It should be noted that at this point there is a signifi- 
cant disagreement between the measurements presented 
here and the data published by the NA57 collaboration 
[a,[3J|. Even though the NA57 data follow the same trend 
in the energy dependence, they are systematically higher 
than the NA49 results 



64j . This discrepancy is generally 



of the order of 1 - 2.5 standard deviations with the only 
exception of the S + measurements at A0A GeV. The 
measured particle ratios, on the other hand, show a good 
agreement between the two experiments. Despite inten- 
sive discussions between both collaborations, the origin 
of the discrepancies is not yet found. 

The energy dependence of the antibaryon/baryon ra- 
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FIG. 13: The total multiplicities of A (a), A (b), S" (c), 
and H + (d) in central Pb+Pb and Au+Au collisions as a 
function of Js BN . The systematic errors are represented 
by the gray boxes. Also shown are AGS data [3, l3a . |45| |. 
as well as calculations with string hadronic models (HSD, 
UrQMD 1.3 04], H, OH) and with a statistical hadron gas 
model (SHM(A) [U). 



and 158,4 GeV. The difference between UrQMD1.3 and 
the statistical model is more pronounced for S~ and S + 
(see Fig. [T3"b and d). While the data points at SPS en- 
ergies are above the UrQMD 1.3 calculation by a factor 
of rs 2, the statistical model fit provides a qualitative de- 
scription of the measurement, although the agreement is 
not perfect. 
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tios _R(B/B), measured at mid-rapidity, are compared for 
protons, A, S, and il in Fig.|T2l The ratios exhibit a rapid 
rise for all particle species over several orders of magni- 
tude in the SPS energy range and converge towards val- 
ues close to 1 at RHIC energies. There is a distinct hier- 
archy of the ratios, depending on the strangeness content 
of the baryons: 

R(n+/n-) > R(Z+/~-) > R(A/A) > R(p/p). 

Also, the energy dependence in the SPS region gets 
slightly weaker with increasing strangeness. The B/B ra- 
tios at mid-rapidity directly reflect the drastic change in 
the net-baryon number. However, the sensitivity depends 
to some extent on the valence quark content of the baryon 
which is thus responsible for the observed hierarchy. 

The total multiplicities, as determined from the 
dN/dy spectra shown in Fig. [5J are compiled in Fig. Q2] 
together with AGS data where available [9|, (3f| |45j . The 
total multiplicities of A and S _ increase quite rapidly 
at lower energies, while from ^s NN sa 8 GeV on they 
rise only moderately with energy. A and S + , on the 
other hand, exhibit a continuous fast increase with beam 
energy. The measurements are confronted with several 
hadronic models. In Fig.[i"3"b calculations with the string 
hadronic models HSD [53J] and UrQMD1.3 [H[ for (A) 
as a function of y/s NN are shown, as well as results from 
a fit with a statistical hadron gas model [29| (SHM(A)). 
All three models are able to describe the data satisfacto- 
rily. A similar picture is observed for (A) (see Fig. IT3"b). 
although the fit with the statistical hadron gas model 
seems to overpredict the measurements at i^boam = 80A 



FIG. 14: The total multiplicities of A (a), A (b), H" (c), 
and H + (d) divided by the total pion multiplicities ((7r) = 
1.5 ((tt + ) + (7i~~))) in central Pb+Pb and Au+Au collisions as 
a function of ^/s NN . The systematic errors are re p resente d by 
the gray boxes. Also shown are AGS data [3. I35J l45l. [54. [5a] . 
as well as calculations with string hadronic models (HSD, 
UrQMD 1.3 [4ll [43 . [53 ]') and a statistical hadron gas model 
(SHM(A) ^). 

In Fig. Q3] the total multiplicities of hyperons divided 
by the total number of pions (w) = 1.5 ((tt + ) + (7r~)) 
are compiled. The (A)/(ir) and the (S~)/(7r) ratios have 
distinct maxima in the region ^/s NN =7-9 GeV, while 
the (A) /(it) and (E + )/(ir) ratios increase monotonously 
with energy. The comparison to the string hadronic 
model results of HSD and of UrQMD1.3 reveals a sig- 
nificant disagreement with the (A)/{ir) and (A)/(tt) ra- 
tios (see Fig. [T4"k and b), which is not present in the 
A and A multiplicities alone as shown in Fig. [T3k and 
b. This is a reflection of the fact that these models 
overpredict the pi qn p roduction at top AGS and lower 
SPS energies [42|, [62| . Hence, the disagreement with 
the (S-)/(tt) and (S + )/(tt) ratios (Fig. [Hfc and d) is 
even more pronounced than for the S multiplicities alone. 
The statistical hadron gas model approach provides over- 
all a better description of the measured particle ratios 
than UrQMD 1.3. However, the {A)/(tt) ratio is clearly 
overestimated at higher energies by SHM(A), while the 
fit results from this model are slightly below the data 



points for (s )/(tt) and (E + }/(n) for 



< 17.3 GeV. 



In [63j it was argued that a statistical model approach 
predicts different positions of the maxima in the energy 
dependence of (A)/(ir) (ys^fmax) = 5.1 GeV) and of 
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FIG. 15: The rapidity densities AN /Ay at mid-rapidity of A 
(a), A (b), H~ (c), and H + (d) divided by the pion rapidity 
densities (tt — 1.5 (% + + 7r - )) in central Pb+Pb and Au+Au 
collisions as a function of ^/s NN . The systematic errors are 
represented by the gray boxes. Also shown are NA57 [al34l . 
AGS [H, [H, 1H, Hl| , and RHIC H3, H H HI [H, OaTHl 
data, as well as calculations with string hadronic models 
(HSD, UrQMD1.3_E|, 5J (H) and a statistical hadron gas 
model (SHM(B) [1^)7 



(S _ )/(7r~) (^/s NN (max) = 10.2 GeV). However, the ex- 
isting measurements do not allow to determine the exact 
positions of the maxima with the required precision in 
order to establish a significant difference. For this pur- 
pose also more data at lower energies (^/s NN < 6 GeV) 
with high precision would be required. 

Qualitatively the same picture emerges when the ra- 
tios of the mid-rapidity yields are studied instead of the 
ratios of total yields, as shown in Fig.[T51together with re- 
sults from RHIC experiments. Again, the string hadronic 
models HSD and UrQMD1.3 fail to match the S~ /tt and 
S+/-7T ratios, even though a reasonable description of the 
A/ir and A/7T ratios at SPS energies is achieved, and sta- 
tistical models provide generally a better description. As 
an alternative implementation of the statistical hadron 
gas model here the one by 0, 11] (SHM(B)) is used. 
While in SHM(A) [29|] a separate fit at each energy to all 
available particle multiplicities is performed by varying 
chemical freeze-out temperature T c h and baryonic chem- 
ical potential /i B , the input parameters T c h and /i B in 
SHM(B) [59( are taken from a smooth parametrization 
of the y/s NN -dependence of the original fit results. In ad- 
dition, the model SHM(A) includes a parameter to allow 
for strangeness undersaturation 7 S , which is not present 
in model SHM(B) (i.e. 7 S = 1). Both models use a grand 
canonical ensemble for the results shown here. In the case 
of model SHM(B) an additional correction by a canonical 
suppression factor is applied. However, for central A+A 
collisions this correction is only effective at AGS energies 
(Js^ < 5 GeV) [HI. Thus, SHM(B) provides a baseline 
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FIG. 16: The total multiplicities of A (a) and A (b) divided 
by the total pion multiplicities ({tv} = 1.5 {(tt + } + (i" - ))) for 
central Pb+Pb and Au+Au collisions as a function of ,/s NN . 



The systematic errors are represented by the gray boxes. Also 
shown are AGS data [y, 13a . l4a . 154 . |55J . measurements for 
p+p collisions by other experiments (open squares) [6fj|. [61|. 
as well as a calulation with the UrQMD1.3 model (dashed 
line: Pb+Pb, dotted line p+p). The solid line represents a 
parametrization of the p+p data (see text). The enhance- 
ments relative to the p+p parametrization are shown in pan- 
els (c) and (d). The gray boxes denote the uncertainty of the 
p+p reference parametrization. The dashed line represents 
the values from the UrQMD1.3 model. 



defining the state of maximal chemical equilibrium that is 
attainable. However, the parametrization that provides 
the basis of SHM(B) has been tuned to fit mid-rapidity 
ratios, while the fits with SHM(A) have been done for 
total multiplicities which complicates a direct compari- 
son between the two approaches. SHM(B), as shown in 
Fig. \W[ generally overpredicts all measured mid-rapidity 
ratios at the higher SPS energies (y/s^ = 12 — 17 GeV), 
while at lower SPS and at RHlG energies a satisfac- 
tory agreement is achieved. Therefore, in the data a 
sharper maximum in the energy dependence of the A/7T 
and 'Er /tt ratios is observed than in the model. The 
NA57 results[65J| exibit a similarly shaped energy depen- 
dence. However, the ratios are generally higher than the 
NA49 results. 

The observed maxima in the A/7T and S~/7r ratios oc- 
cur in the same energy range as the observed distinct 
peak in the K + /ir + ratio [2f|. Since the latter can be in- 
terpreted as a signature for the onset of deconfinement, 
the question appears whether the maxima in the A/7T 
and S~/7r ratios can be attributed to the same effect. 
In contrast to the K + , which carry together with the 
K° the bulk of the anti-strange quarks and are thus a 
relatively direct measure of the strangeness production, 
the interpretation of the strange baryons is complicated 
by the fact that their sensitivity to the strangeness pro- 
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duction is strongly modified by the energy dependent 
baryon number distributions. At low energies, with high 
baryonic chemical potential, the production of baryons is 
favoured and more strange quarks will end up in A and 
S _ , compared to higher energies where strange quarks 
might predominantly be contained in K~ and K . This 
is underlined by the fact that the statistical model ap- 
proaches, which reflect the dependence of particle yields 
on fi B , provide a relatively good describtion of the data. 
Whether the remaining discrepancies between SHM(B) 
and the mid-rapidity ratios at 80A and 158A GeV (see 
Fig. [T5j panels (a) and (c)) might be attributed to the 
onset of deconfincment can in the light of the systematic 
uncertainties not be definitely answered. However, one 
should keep in mind that the mid-rapidity A/n and E+ /w 
ratios are also strongly affected by the rapid change of 
the shape of the A and S - rapidity distributions with 
energy. This effect will cause a more pronounced energy 
dependence of the mid-rapidity ratios in comparison to 
the 4-7T ratios, which in principle cannot be described by 
statistical models. 

The (A)/(7r) and (A)/(tt) ratios, as measured in cen- 
tral nucleus-nucleus collisions, are compared to data ob- 
tained in p+p collisions [60, [61| in Fig. [THI The p+p 
measurements were parametrized by a fit function. For 
(A)/(7r) the following function with the fit parameters a, 
b, and c was used: 

(A)/(tt)(p+ p ) = c[l-exp(-(VS— s/io)/a)+6(Vi- 



~o)] 

Here, y/so denotes the threshold center-of-mass energy. 
The result of the fit is displayed in Fig. [T6k . It pro- 
vides a reasonable description of the available data in 
the energy range of ^/s^ < 20 GeV. Similarly, the en- 
ergy dependence of the (A)/(ir) ratio was parametrized 
by a straight line. However, the existing measurements 
are much less precise than in the (A)/(ir) case. Based 
on these parametrizations, the energy dependence of an 
enhancement factor E relative to p+p, defined as 



E=(N)/(7r)\ A+A /(N)/(w)\ p+p 



(5) 



can be determined. As shown in Fig. 116b . the enhance- 
ment factor for A exhibits a clear increase from a factor of 
«2 to >3 towards lower energies. For ^/s NN < 4 GeV the 
AGS measurement of [9( suggests an even more dramatic 
rise towards very low energies. For A the enhancement 
is of the order of rs 2, without any significant energy 
dependence in the range covered by the data. While the 
UrQMD1.3 model qualitatively reproduces the energy de- 
pendence of the A-enhancement, it fails to describe the 
enhancement of A. In fact, the model rather predicts a 
A-suppression, which is mainly due to the fact that the 
(A)/(7r) ratio in p+p reactions is grossly overestimated 
(see dotted line in Fig. HHb). Since the net-baryon den- 
sity is largest around ^/s NN = 5 GeV, the production 
of strange baryons exhibits a pronounced maximum at 
these energies. This effect is described by all hadronic 
models considered here and consequently the A/-7r-ratios 



are well reproduced (Figs. Q3] and [T5|) . Moreover, the 
energy dependence of the A-enhancement seems to be af- 
fected by the redistribution of the baryon number, which 
is suggested by the fact that UrQMD1.3 gives a similar 
increase towards low energies. In comparison, the dou- 
bly strange !B _ is less sensitive to the baryon number 
density and more to the overall strangeness production, 
which may explain why string hadronic models fail to de- 
scribe the data. For the corresponding antiparticles this 
argument applies even more strongly. Whether the 3- 
enhancement also increases towards low energies, similar 
to the A, can currently not be decided due to the lack of 
precise reference data in p+p at lower energies. 



VI. SUMMARY 

A systematic study of the energy dependence of A, A, 
E — , and S + production in central Pb+Pb reactions at 
SPS energies is presented. 

The shape of the rrit spectra exhibits only a weak de- 
pendence on beam energy, which is also reflected in the 
moderate increase of (mt)-mo towards the higher RHIC 
energies. A similar behavior was also observed for pions, 
kaons, and protons. For these particles a sudden change 
in the energy dependence around ^/s NN =7-8 GeV was 
found in addition. Due to the lack of data at lower ener- 
gies it currently cannot be established whether a similar 
feature is present in the energy dependence of (tnt) — too 
for hyperons. There is an indication for a slightly weaker 
energy dependence of (mt) — too f° r A than for A, the 
values for A being above the ones for A. Generally, the 
measured (rat) — too is higher for all investigated particle 
species than what is predicted by the UrQMD1.3 model. 

For A and S~ rapidity spectra a clear change of the 
shape is observed. The almost Gaussian like spectral 
form develops a plateau around mid-rapidity towards 
higher energies, reflecting the change of the longitudi- 
nal distribution of the net-baryon number. The rapidity 
spectra of A and S + , on the other hand, can be described 
by single Gaussians at all investigated energies, whose a 
increases monotonically with energy. 

Also the energy dependence of the total yields shows 
a distinct difference between baryons and anti-baryons. 
While for the A and S + multiplicities a continuous rapid 
rise with beam energy is observed, the increase of the 
A and E~ yields is clearly weaker above y/s^ = 7 - 
8 GeV than below. This difference gets even more pro- 
nounced when dividing the total multiplicities of the hy- 
perons by those of pions. The energy dependence of the 
(A)/(ir) and (EE+)/(7r) ratios exhibits significant maxima 
in the region 5 < ^/s NN < 10 GeV, while the (A)/(tt) 
and (S + )/(7r) ratios increase monotonically. The total 
multiplicities of A and A are well described by the string 
hadronic UrQMD1.3 model. However, S~ and S + mul- 
tiplicities are underpredicted by factors of 2 - 3 at SPS 
energies. A better overall description of all measured 
yields is provided by statistical hadron gas models. 
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